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Abstract 

The basic concepts and hypotheses of Newtonian Cosmology necessary for a 
consistent treatment of the averaged cosmological dynamics are formulated and 
discussed in details. The space-time, space, time and ensemble averages for the cos- 
mological fluid fields are defined and analyzed with a special attention paid to their 
analytic properties. It is shown that all averaging procedures require an arrange- 
ment for a standard measurement device with the same measurement time interval 
and the same space region determined by the measurement device resolution to be 
prescribed to each position and each moment of time throughout a cosmological 
fluid configuration. The formulae for averaging out the partial derivatives of the av- 
eraged cosmological fluid fields and the main formula for averaging out the material 
derivatives have been proved. The full system of the averaged Navier-Stokes-Poisson 
equations in terms of the fluid kinematic quantities is derived. 



1 Introduction 

In order to describe the dynamics of inhomogeneous and anisotropic Newtonian universes 
one needs to establish a set of averaged Navier-Stokes-Poisson equations in kinematic 
quantities. To approach this problem the basic concepts and hypotheses of Newtonian 
Cosmology necessary for a consistent treatment of averaged cosmological dynamics are 
formulated and discussed in details. The space-time, space, time and ensemble averages 
for the cosmological fluid fields are defined and analyzed with a special attention paid for 
their analytic properties. It is shown that all averaging procedures require an arrangement 
for a standard measurement device with the same measurement time interval and the same 
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space region determined by the measurement device resolution to be prescribed to each 
position and each moment of time throughout a cosmological fluid configuration. The 
formulae for averaging out the partial derivatives of the average fluid flelds and the main 
formula for averaging out the material derivatives have been proved. The full set of the 
averaged Navier-Stokes-Poisson equations in terms of the fluid kinematic quantities is 
derived. 

The structure of the paper is as follows. The Averaging problem in general relativity 
which has originated in cosmology when it has been realized that a modiflcation of the 
Einstein equations by smoothing them out in some appropriate sense should be estab- 
lished for a description of the large-scale structure of the Universe is discussed in Chapter 
2. It is shown that Newtonian cosmology also has the Averaging problem, since due to 
the nonhnearity of the Navier-Stokes-Poisson equations one needs to prove that the av- 
eraged Navier-Stokes-Poisson equations have a solution for a homogeneous and isotropic 
Newtonian universe which can be obtained upon averaging out a locally inhomogcncous 
and/or anisotropic cosmological fluid configuration representing evolving matter struc- 
tures. Chapter 3 shows that an averaged description of a Newtonian universe means, 
from the physical point of view, that the dynamics of the Newtonian universe is written 
directly in terms of observationally relevant measurable quantities. The Newtonian cos- 
mological macroscopic hypothesis which assumes the existence of two essentially distinct 
types of temporal and spatial scales, namely, the local scales of cosmological fiuid fiuctu- 
ations and cosmological mean scales, is formulated and discussed in details in Chapter 4. 
A classification of Newtonian cosmologies is given, in particular, on the basis of this hy- 
pothesis. Chapter 5 gives the definitions of the space-time averages of cosmological fiuid 
fields and their correlations and discusses how one can construct an average fluid field by 
arranging a covering of the Newtonian space-time by a standard measurement device at 
every moment of time and at every position throughout a cosmological fluid conflguration. 
The Reynolds conditions necessary for any kinds of averages arc formulated and proved 
for the space-time averages in Chapter 6. Properties and conditions on the space-time 
correlation functions are proved in Chapter 7. Chapter 8 is devoted to the deflnition 
and analysis of time averages and time correlation functions. It is shown, in particular, 
that time averages over an inflnite time interval docs not depend on time. Chapter 9 is 
devoted to the deflnition and analysis of space averages and space correlation functions. 
It is shown that space averages over the whole space are constant with respect to the 
space position. The formula for the space-time averaging out of the material derivatives 
is proved in Chapter 10 and the form of the space-time averaged Navier- Stokes operator is 
derived. In Chapter 11a possibility to use the Reynolds transport theorem for a deflnition 
of space averages is analyzed and it is shown that one cannot take this approach for the 
proper treatment of space averages. In Chapter 12 the ensemble averages and ensemble 
correlation fimctions arc deflned and the statistical treatment of cosmological fluid flelds is 
discussed. The Hypothesis of the statistical nature of Newtonian universes is also formu- 
lated here. Properties and conditions on the ensemble averages and ensemble correlation 
functions, such as the Reynolds conditions, the formula for the ensemble averaging out of 
the material derivatives and the form the ensemble averaged Navier- Stokes operator, are 
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proved in Chapter 13. The Ergodic hypothesis of Newtonian cosmology is discussed and 
formulated in Chapter 14. The full set of the averaged Navier-Stokes-Poisson equations 
in terms of the fluid kinematic quantities valid for the space-time and ensemble averages 
is derived in Chapter 15. It is pointed out that the system must be supplemented by a 
system of the equations for the single-point second order moments of fluid fields present 
in the averaged Navier-Stokes-Poisson equations. 

The formulae and Sections from the papers I and II [|l[] , will be referred to as (I-XX) 
and I-X, and (II-XX) and II-X, correspondingly. 

Conventions and notations are as follows. All functions / = are defined on 

3-dimensional Euclidean space in the Cartesian coordinates {x*} with Latin space 
indices i,j,k,... running from 1 to 3, and t is the time variable. The Levi-Civita sym- 
bol Eijk is defined as £123 = +1 and e^^^ = +1 and 5*-^, 5* and 6ij are the Kronecker 
symbols. The symmetrization of indices of a tensor Tj^i, = T'^-j^{x\t) is denoted by 
round brackets, T^^^^ = ^ (Tj^ + T^^), and the antisymmetrization by square brackets, 
T^k] ~ \ i'^jk ~ "^kj)- ^ partial derivative of Tj^ with respect to a spatial coordinate x* or 
time t is denoted either by comma, or by the standard calculus notation, Tj^^^ = dT^/^/dt 
and Tjf^ i = dTjj^/dx^ , the fluid velocity is = u'^{x^ ,t), and the material (total) derivative 
is T]^ = dT^f^/dt = dT^i^/dt + u^dTj^/dx^ = T];.^ + m'T];.^. Average values are shown by 

either the angular brackets, or by a bar over a tensor, (Tj^) = with the corresponding 
subscripts to distinguish between space-time, time, space and ensemble averages. The 
Newton gravitational constant, the velocity of light and the Newtonian cosmological con- 
stant are G, c and A, respectively. 

2 The Averaging Problem in Newtonian Cosmology 

Our Universe on the largest scales is assumed to be homogeneous and isotropic though 
locally, at the scales of stars, galaxies and even clusters of galaxies, it is inhomogeneous 
because of the matter condensation into these discrete structures. Observational data def- 
initely shows a high degree of isotropy which is supported specifically by measurements of 
the temperature of the cosmic background radiation. The large-scale spatial homogeneity 
of the Universe up to scale of order 100-300 Mpc, however, still remains an issue which 
cannot be considered as undoubtedly proved on the basis of cosmological observations. 
In addition to uncertainties in the measured data, difficulties in proving this fundamental 
hypothesis originate also in the way how this data are treated. In particular, some theo- 
retical models for analysis of observational data are based on an assumption about such 
a homogeneity Therefore, the general relativistic Friedmann-Lemaitre- Robertson- 
Walker (FLRW) cosmological models which are homogeneous and isotropic have 
not been proved as yet to be precisely realistic to represent the large-scale structure of our 
Universe even for the present epoch. The same is valid for the Newtonian cosmological 
models, see Section II-9. 

The observations of the large-scale characteristics of cosmological matter and radiation 
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do measure their averaged values over space regions and time intervals. From physical 
point of view, one can think of the large-scale cosmological models as representing an 
averaged picture of the Universe which is highly inhomogeneous on the smaller scales 
characteristic of typical matter structures. From the mathematical point of view one can 
think of a cosmological model which is inhomogeneous on some scales to mimic local 
gravitating matter condensations and which, at the same time, is homogeneous on a large 
scale in some averaged sense, that is, after application of a smoothing procedure to the 
physical fields and their governing equations. Such a cosmological model would pretend 
to be more adequate to describe the local dynamics of the real locally inhomogeneous Uni- 
verse and its large-scale homogeneous and/or isotropic dynamics. Then the Newtonian 
cosmological principle, see Section II-6, is related to the large-scale space-time structure 
of a locally inhomogeneous universe averaged out over the space regions of a size much 
larger than a typical length of inhomogeneities, say, the curvature radius of local matter 
structures, and much less than a typical large-scale cosmological length, say, the curva- 
ture radius of the Universe due to the mean matter density, and the time intervals much 
larger than a typical time of change of local cosmological variables governing the dynam- 
ics of the matter structures and much less than a typical time of change of cosmological 
variables governing the evolution of the Universe as a whole. The Averaging problem in 
the general relativistic cosmology calls for a development of such a theoretical physically 
reasonable and mathematically rigorous framework where one can formulate a realistic 
inhomogeneous cosmological model which is homogeneous and/or isotropic on average. 
The first attempt to carry out a space-time volume averaging out the Einstein equations 
had been made by Shirokov and Fisher . The Averaging problem has been first appar- 
ently formulated by Ellis |]T0| in the framework of general relativistic cosmology (see |1T 



for a comprehensive review and references). The first covariant approach for averaging 
out the Einstein equations over space-time regions has been developed by Zalaletdinov 



12 1- 1 14 1 in the framework of macroscopic gravity, a classical macroscopic theory of grav- 



ity with the space-time volume averaged (macroscopic) Einstein equations together with 
equations for the gravitational correlation functions governing the dynamics of continuous 
macroscopic matter distributions. The Einstein equations are considered as the micro- 
scopic field equations suitable for description of point-like matter distributions [0, |T5| . 
The gravitational correlations have been shown to be very important for the dynamics 
of the averaged (macroscopic) gravitational field, in particular, when all correlations are 
assumed to vanish the macroscopic Einstein equations become the Einstein equations of 
general relativity. 

Thus, according to the Averaging problem a modification of the Einstein equations by 
smoothing them out in some appropriate sense should be established for a description of 
the large-scale structure of the Universe and an adequate large-scale macroscopic cosmo- 
logical model should be found as a solution to these equations. The Einstein equations 
of general relativity then have to describe the local dynamics of the Universe by a suit- 
able cosmological model. Upon an appropriate averaging of this cosmological model one 
must obtain the macroscopic cosmological model which is homogeneous and/or isotopic 
on some large cosmological scale. 
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The same problem exists for Newtonian cosmology. A distribution of stars, galaxies 
and clusters of galaxies in the real Universe can be represented by a locally inhomoge- 
neous distribution of the cosmological fluid which is governed either by the Navier-Stokes- 
Poisson equations (II-4), (II-8), (II-9) and (11-10), or the system of Navier-Stokes-Poisson 
equations in terms of kinematic quantities (II-4), (II-8), (11-12), (II-51)-(II-56), (11-59)- 
(11-61), (I-64)-(II-66), (11-69), (11-28) and (11-70). Assuming the Newtonian cosmological 
principle, see Section II-6, valid for the large-scale structure of a Newtonian universe, one 
questions now whether or not the suitably averaged Navier-Stokes-Poisson equations or 
the Navier-Stokes-Poisson equations in terms of kinematic quantities have a solution for 
an homogeneous and isotropic Newtonian universe such as this large-scale homogeneous 
and isotropic cosmological fluid configuration can be obtained upon averaging out the lo- 
cally inhomogeneous and / or anisotropic distribution of a cosmological fluid configuration 
representing evolving matter structures. 

The Newtonian averaging problem Given a locally inhomogeneous and/ or anisotro- 
pic distribution of a cosmological fluid configuration representing the evolving cosmologi- 
cal matter structures as a solution to the Navier-Stokes-Poisson equations (11-4), (^^-8), 
(11-9) and (11-10), if there exists a homogeneous and/or isotropic Newtonian universe as 
a solution to the suitably averaged Navier-Stokes-Poisson equations such as the homoge- 
neous and isotropic cosmological fluid configuration can be obtained upon averaging out 
this locally inhomogeneous and/or anisotropic cosmological fluid configuration. Then this 
locally inhomogeneous and/ or anisotropic Newtonian universe will satisfy the Newtonian 
cosmological principle for its large-scale structure. 

The importance of study of this difficult challenging problem is evident for both New- 
tonian and general relativistic cosmologies. At the very fundamental level, the Averaging 
problem questions the validity of the Newtonian cosmological principle, see Section II-6, 
as a selection principle for a realistic locally inhomogeneous and/or anisotropic model of 
the Universe to be globally homogeneous and isotropic on its largest scales. From the 
mathematical point of view, in general relativity due to the nonlinearity of the Einstein 
equations one cannot expect this is to be true for an arbitrary inhomogeneous cosmologi- 
cal model. Taking an average of the Einstein field operator for the metric tensor serving as 
the gravitational potential in general relativity does not commute with taking an average 
of the metric tensor as it has been first pointed out by Shirokov and Fisher [0] (see also 



rO[-|T2[). In fiuid mechanics the same argument is valid for the Navier-Stokes-Poisson 
equations as it has been known after Reynolds |T^. Reynolds was the first to realize that 
the nonlinearity of the Navier-Stokes equation (1-84) actually contains a hidden key inside 
it to the understanding of fiuid motion through definition of the correlation functions of 
fluid velocity and density for analytical study of the dynamics of evolving fluids by means 
of a suitable time and/or space averaging procedure . In Newtonian cosmology this im- 
portant argument has never been formulated explicitly. Indeed, taking an average value 
of the Navier-Stokes field operator, the left-hand side of Eq. (11-10), does not produce 
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the material derivative of an average of the fluid velocity, 




d{u') 



dt 




d{u') 



(1) 



where (■) stands for an averaging operator. Thus, similar to the Averaging problem in 
general relativity one cannot expect that an arbitrary inhomogeneous Newtonian universe 
would provide a solution to the Newtonian averaging problem. A thorough analysis is 
required to understand if there exist some locally inhomogeneous Newtonian universes 
satisfying the Newtonian cosmological principle on their largest scales. 

3 Averaging and Space-time Measurements 

For an inhomogeneous and/or anisotropic Newtonian universe, the spatial and temporal 
dependence of the instant values of the cosmological fluid velocity = u'^{x^,t) are very 
complex. Moreover, if the cosmological fluid flow is turbulent the values of the fluid 
velocity fleld are known to be different each time even if the initial conditions have been 
set the same for each fluid conflguration. Due to this extreme uncertainty and sharp 
changes in space and time of the cosmological fluid velocity fleld the density 

fleld p{x^,t) and the pressure fleld p{x^,t), it is necessary to use an averaging procedure 
to introduce more regular mean values of the fluid characteristic quantities instead of real 
fast and sharp varying dynamical fluid flelds. These averaged flelds may then be studied 
analytically by means of usual methods of mathematical physics. 

Another justiflcation for using the averaged quantities for characterization of the fluid 
motion comes from the measurement theory. Indeed, every real act of measurement of 
the velocity, the density or the pressure of a fluid particle is carried out by an observer for 
some interval of time T = t2 — ti during which the particle is moving to another position 
{^2} — {^1 + ^x^}. The measurement time interval, and, correspondingly, the loca- 
tion where a fluid particle moves during the measurement, depends on the measurement 
device and conditions. To reduce a measurement error in the determination of a fluid 
fleld below an acceptable or required value, a measurement may take a considerable time 
and the influence of possible random fluctuations in the measured fluid quantity value 
is then reduced as necessary. The measured value of the fluid fleld is an average over a 
measurement time interval T and over a space region S G which is determined by 
the precision of the position resolution of the measurement device. All fluctuations in the 
quantity under measurement with the temporal variations much less than the measure- 
ment time and with the spatial variations of much smaller scales than the measurement 
device resolution scale are completely suppressed without any trace of them in the mea- 
sured average valueQ. Therefore, a formulation of Newtonian cosmology in terms of space 

^ One should bear in mind, however, that it only concerns the linear fluctuations of fluid charateristics. 
Average values of quadratic and higher-order combinations of fluctuating fields do not vanish in general 
to bring the fluid moment, or correlation, functions which are essential for understanding of the dynamics 
of moving fluids , ||l8| , see Section |[ 
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and time averages of the cosmological fluid velocity, density and pressure means, from the 
physical point of view, that the dynamics of the Newtonian universe is written directly in 
terms of observationally relevant variables, and all possible physical effects such as fluid 
correlations have been introduced into the theoretical framework during derivation of such 
an averaged theory. 



4 Averaging out a Cosmological Fluid Field 

A physical picture of an inhomogeneous Newtonian universe is that any fluid particle 
has different values of its instant velocity t), the density p{x\t) and the pressure 

t) as it moves moving along its path. All other fluid particles have different values of 
the instant velocity, the density and the pressure as they move along their paths and as 
each compared with others. In order to describe now a cosmological fluid configuration 
in terms of the mean values of fluid characteristics one should be able to construct the 
average fields for the fluid velocity the fluid density (p(x^,t)) and the fluid 

pressure (p(x'^,t)). A cosmological fluid configuration with the average characteristics 
may have a different dynamics as compared with the original fluid configuration as far as 
the average fields are not equal to the original fluid fields, 

{u\x^,t))^u\x^,t), {p{x\t)) ^ p{x\t), {p{x\t))^p{x\t). (2) 

An equation of motion = x^{C,\t) for a fluid particle moving in an inhomogeneous 
self-gravitating fluid configuration and an equation of motion X* = , t) for the same 
fluid particle moving in an averaged self-gravitating fluid configuration are expected to 
be different because of different solutions of the initial value problems for the the fluid 
velocity field u\x^,t), see (I-IO) and (11-31), 

= u%x^,t), x\0)=C, (3) 

at 

and for the average fluid velocity field {u\x\t)) (Q) 

{u\x\t)), X\0) = e, (4) 



dt 



even if the fluid particle has been given the same initial values x^{0) = and X*(0) = 
in both fluid configurations. So the paths of the same fluid particle as it is moving within 
an inhomogeneous cosmological fluid distribution and within the corresponding averaged 
cosmological fluid distribution have to be different in general. 

The question of the definition of the average values of classical physical fields is a 
delicate one and it has a long history in the classical theory of fields (for a discussion and 
references see |17|, |^0| for fiuid mechanics, |21]- for the classical macroscopic 



electrodynamics, and |TT||, |jT^, p5[ for general relativity). In practice, as it has been 



pointed out above one usually deals with a space-time average values as resulting from 
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observations and measurements. Let us consider a quantity / which characterizes the 
fluid in the framework of the field description. Generally, for an inhomogeneous and/or 
anisotropic cosmological fiuid configuration evolving in time such a quantity / is a function 
of the Eulerian variables / = with nonvanishing temporal and spatial derivatives 

which define a temporally and spatially varying field f{x^,t) for every fiuid particle. To 
provide a nontrivial averaged cosmological dynamics the cosmological fiuid field f{x\t) 
is assumed to satisfy the following physical hypothesis: 

The Newtonian cosmological macroscopic hypothesis Given a measurement time 
interval T and a space region S G with a volume Vs determined by the measurement 
device resolution, a cosmological fluid field f{x\t) is supposed to have, at least, either 
two essentially different temporal variation scales (At, Lt) , or two essentially different 
spatial variation scales {Xs,Ls), or both pairs of scales simultaneously, such as 

At « T « Lt, Xs « V^^^^ « Ls (6) 

where At and A5 are the temporal and spatial scales of local fluid field fluctuations, and 
Lrp and Ls are the temporal and spatial mean cosmological scales, correspondingly. 

Since there is in fact a hierarchy of distinct physical scales in our Universe |jlO|, the 



validity of the Newtonian cosmological macroscopic hypothesis for a fiuid field 
should be checked for each pair of scales of interest. 

Hypotheses analogous to the Newtonian cosmological macroscopic hypothesis (^) are 
always assumed, very often implicitly, for every time, space, or space-time averaging 
scheme (see [|T7|, [|T^ for fiuid mechanics, ||2T|-[^ for the classical macroscopic electro- 



dynamics, and [0, p5| for general relativity). The physical picture of a Newtonian 
universe satisfying the Newtonian cosmological macroscopic hypothesis (^ corresponds 
to our intuitive understanding of our Universe based on observations and space exper- 
iments. The Universe is seen now as highly inhomogeneous on the scales of different 
matter structures which are represented by temporal and/or spatial scales of local fluid 
fleld fluctuations At and A5. On its largest scales the Universe shows a deflnite tendency 
to homogeneity and isotropy which are represented by temporal and/or spatial mean cos- 
mological scales, Lt and Ls- On some stages of the evolution of our Universe one can 
expect, however, that the hypothesis (|^) may not be valid or has to be modifled. For 
instance, a consideration of the Universe evolution at early times when it might have 
had a turbulent regime, the cosmological fluid may have more complicated fluctuation 
frequency proflles. In this situation, using a mean fleld picture and the corresponding 
averages should thoroughly investigated to verify whether or not such a consideration is 
physically relevant and adequate, see Sections II-6 and ^ 

In dependence of whether the Newtonian cosmological macroscopic hypothesis (^ 
valid or not for a cosmological fluid distribution, one can deflne three possible classes of 
Newtonian universes. 
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Class LC (Locally inhomogeneous and/or anisotropic Newtonian universes 
with the large-scale cosmological dynamics) Such Newtonian universes satisfy the 
Newtonian cosmological macroscopic hypothesis ^ and the averaged cosmological fluid 
fields (/(x*, t)) has a nontrivial dynamics with either the characteristic temporal variation 
scale Xrp, or the spatial variation scale A5, or with both scales of local fluid field fluctuations 
with the condition (|^ being held. 

This class of Newtonian universes, if it is not empty, would provide a solution to 
the Averaging problem, see Section |^. One of the Newtonian cosmological models of this 
class may be a realistic cosmological model of our Universe in the framework of Newtonian 
cosmology. 

Class L (Locally inhomogeneous and/or anisotropic Newtonian universes with- 
out the large-scale cosmological dynamics) // a cosmological fluid configuration 
does not possess any temporal and/or spatial mean cosmological scales, Lt and Ls, when 
the averages of the cosmological fluid fields vanish and there is no large-scale cosmological 
regime of the evolution of such Newtonian universes, that is, 

{u\x\t)) = 0, {p{x\t)) = 0, {p{x\t))=0. (7) 

This class of Newtonian universes has clearly no physical interest when one is look- 
ing for a realistic cosmological model. However, it is of interest in connection with the 
Averaging problem since, if it is found, one may be able to formulate a set of definite 
constraints on the local matter distribution necessary for a Newtonian universe to belong 
to the Class LC. 

Class C (Newtonian universes with the large-scale cosmological dynamics 
only) // a cosmological fluid configuration does not possess any temporal and/or spa- 
tial scales of local fluid field fluctuations Xt and Xs , when no nontrivial averaged regime 
is possible and the average cosmological fluid fields are equal to the original fluid fields, 
that is, 

{u\x\t)) = u\x\t), {p{x\t)) = p{x\t), {p{x\t))=p{x\t). (8) 

This class of Newtonian universes does not seem to be realistic since our Universe does 
have temporal and/or spatial scales of local fluid field fluctuations At and As, though 
a locally inhomogeneous matter distribution. Again, it is of interest in connection with 
the Averaging problem since, if this Class is found, one may be able formulate a set of 
definite constraints on the large-scale cosmological dynamics necessary for a Newtonian 
universe to belong to the Class LC. 

As a fundamental physical consequence of the Newtonian cosmological macroscopic 
hypothesis (|^), an averaging region S" x T of a space-time measurement is considered 
effectively as a single "point" for the averaged cosmological fluid field {f{x\t)). Such 



9 



regions have been called "physically infinitesimally small" by Lorentz |^T| and the results 
of the measurements must be insensitive to a choice of a reference point to which 

the average {f{x\t)) is prescribed, t eT, {x''} G S. Under the Newtonian cosmological 
macroscopic hypothesis (^) any measurements of cosmological fluid quantities over a time 
interval T and a space region S (Z with a volume V5 determined by the measurement 
device resolution guarantee that the errors made while performing such measurements 
in different positions {x^} within the region S, {x^} G 5", and different instants of time 

— 1/3 

t within the interval T,tET, are always less than 0{\t/T) and 0{\s/Vg ) (see for 



a discussion and references therein) . The possibility of choosing the averaging time interval 
T and a space region S G with a volume Vs to be intermediate between the fluctuating 
and mean fluid field scales, Ay, A5 and Lt, Ls, assumes that the cosmological fluid motion 
may be resolved into a smoothly varying averaged motion with a very irregular fluctuating 
motion superimposed on it. There is a considerable gap between the frequency range 
characteristic for the mean motion and fluctuating motion. 

In studying the averaged (mean) cosmological fluid configurations there are three 
following types of averages of a cosmological fluid field /(x*, t) of interest: a time average 
{f{x\t))rp, a volume space average {f{x\t))g and an ensemble, or statistical, average 
{f{x\ t)) ^ . The first two types of averages can be naturally defined and treated together 
as a space-time average {f{x\t))gj, . 



5 The Space-time Averages and Correlations of a 
Cosmological Fluid 

Given a measurement time interval T and a space region S G E^ determined by the 
measurement device resolution, a space-time average of the fluid field / = f{x\t) is 



defined fll, , pg, [||] as 

The space-time average of a fluid field / 

where Vs is 3-volume of the space region S, 



fix' + x'\t + t')dt'dV' 



Vs = I dV 

s 



(9) 



(10) 



The space-time average (P) is a well-defined function of a reference moment of the 
absolute time t, t E T, and of a reference position {x''}, {x^} G S, with the space-time 

^ Sometimes a continuous weighing function vanishing outside the space-time region S x T is used 
in the definition (j^) where the integration is then carried out over the whole four-dimensional manifold 
defined for the physical confuguration under consideration. Such a weighting function is completely 
determined by the properties of a measurement device and the experimental conditions. 
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average (/(x*, t))gj, being prescribed to (x*, t). If the Newtonian cosmological macroscopic 
hypothesis (j^) is satisfied, the choice of the reference time and position is arbitrary within 
the averaging region T x S. The average value is generally a functional of T and Vs. 

Because the Navier-Stokes-Poisson equations are nonlinear one faces necessity to deal 
with averaging products of fiuid fields. The space-time average (|) for a product of 
cosmological fiuid fields / = f{x\t) and h = h{x^,t) defines the so-called correlation 
function. 

The two-point second order space-time correlation, or moment, function of 
fluid fields / and h The two-point second order correlation function is the space-time 
average of a product of f{x\t) and h{y^, s) 



{f{x\t)h{y\s))^^ = ^ I I f{x' + x'\t + t')h{y^+x'^,s + t')dt'dV'. (11) 



's Jt Js 



The correlation function (|TT]) is a well-defined function of the reference moments t and s of 
the absolute time, t,s E T, and of the reference positions {x^} and {y''}, {x'^}, {y''} G S, 
with the space-time average {f{x\t)h{y^,s))gj, being prescribed to the pair of Eulerian 
variables {x\t) and {y\s) and symmetric in these variables 

{f{x\t)h{y^,s))^^ = {f{y^,s)h{x\t))^^ (12) 



From the physical point of view the two-point space-time average { ^A\ j corresponds to 
a simultaneous measurement of / and h by an observer with a measurement device with 
time measurement interval T and a measurement device resolution region S G E'^. As far 
as the Newtonian cosmological macroscopic hypothesis @ is satisfied for the cosmological 
fluid configuration under consideration a choice of the reference times t and s and of the 
reference positions {x^} and {y^} is arbitrary within T and S. The average products 
of fluid fields are known |T^, |T^, to be fundamental characteristic functions 



responsible for essentially nonlinear phenomenon in evolving fiuids, such as, for example, 
turbulence and hydrodynamic instability. If some of such functions do not vanish for a 
fluid conflguration, then the corresponding fluid flelds are said to have correlations. 

The fluid field correlations A cosmological fluid configuration has a central two-point 
second order space-time correlation, or moment, function Cgrp{x\t;y\s) , if there are, 
at least, two cosmological fluid fields f{x\t) and h{x\s) such that 

Cfl{x\t-y\s) = {f{x\t)h{y\s))^^ - {f{x\t))^^{h{y^,s))^^ ^ 0, (13) 

4'^(x^^;y^s) = 4'^(y^s;a;^t) (14) 

in some open space region U G for an interval of time At, x^ , y^ E U, t,s E At. 

One can deflne the central multi-point higher-order space-time correlation functions 
Cg^(a;*, t; y\ s; z'^, r) and so on, whenever it is necessary for analysis of the dynamics of a 
Newtonian universe by making deflnitions similar to (|ll|) and (plSj). 
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In order to characterize the space-time averaged fluid configuration one needs to deter- 
mine now an average field {f{x\t))gj, and a correlation function field {f{x\t)h{y^,s))grp 
for each moments of time t and s at all positions {x^} and {y^}. From the physical point 
of view, that means that it is possible, at least in principle, to carry out the measurements 
of the quantities / and h for all instants of time during the evolution of a fiuid configu- 
ration and at its each space point. This requires additional assumptions concerning the 
measurement interval T and the space region S which are usually made only tacitly, or 
they are supposed to be trivial (see a discussion in [^], [^). This is, however, a very 



significant issue for a correct definition of analytical properties for an average field (see 



121 , IM , [g^ for a detailed discussion). 

The First condition (A covering of the Newtonian space-time by averaging 
time intervals and space regions) A measurement time interval T and a space region 
S must be prescribed at every moment of time t and every point point {x'^} in order to 
define an averaged fluid field. 

Given such a covering by the measurement interval T and the space region S, the set of 
the average values calculated at each moment of time t and at each position {x^} will form 
average fields {f{x\ t))gj, and {f{x\ t)h{y^ , s))^^,. As noted above, this set corresponds to 
the average values of / as they are measured for all instants of time during the evolution 
of a cosmological fiuid configuration at its each space point. 

The Second condition (Typical averaging time intervals and space regions) All 

time intervals T and the space regions S are typical in some defined sense. 

They are usually required to be of the same shape and volume, T = const and Vs = 
const, such as 



at each moment of time t and at each position {s'^}. For a Newtonian universe which 
has the Newtonian space-time geometry, see Sections 11-3 and 11-5, one can always satisfy 
both conditions by placing a measurement device at all positions {x^} in at an initial 
moment of time t = to and by the Lie-dragging of the time interval T and the region 
5* along the congruences of the Cartesian coordinate lines x'^ and the absolute time t 
to get a covering of the manifold by the regions S of the same shape and volume and 
the time intervals of the same length T at every point of the Newtonian space-time. 
Another way of arranging such a covering is to associate a measurement device to all 
positions {x''} in at an initial moment of time t = to with fiuid particles located 
at {x'^} at t = to and transport the device to each other position during the evolution 
of a cosmological fiuid configuration by comoving with the corresponding fiuid particles 
as they evolve. Both prescriptions must bring the same covering, that is, associate a 
standard measurement device with the measurement time interval T and the measurement 
resolution space region S with each point {x^} at each instant of time t throughout the 
cosmological fiuid configuration. Three issues here are of particular importance for notice: 
(^4) Though the second procedure involving transportation along cosmological fiuid 
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paths of moving fluid particles seem to be more physical, one should bear in mind that 
this way of construction of a field of the measurement devices relies on the dynamical 
evolution of a cosmological fluid. In many case it does not cause any complications. How- 
ever, in situations when the dynamics of an inhomogeneous Newtonian universe exhibits 
some specific fluid path configurations, like asymptotic attractors, there might be some 
points in the Newtonian space-time where no measurement device can be brought by this 
procedure. The first procedure is essentially kinematic and relies only on the topological 
and differential structure of the underlying Newtonian space-time manifold for an inho- 
mogeneous Newtonian universe. The Eulerian variables are defined globally on the 
Newtonian space-time manifold, therefore by the Lie-dragging of the measurement devices 
along the coordinate lines (x*, t) they can be determined for all times t and positions {x^}. 
The Lie-dragging itself is a physically well-defined and meaningful procedure since it pre- 
serves physical scales along the coordinates lines (x* , t) and does not affect the manifold 
symmetries and structure. It does fit perfectly for the purpose of the construction of a 
covering by averaging regions without any possible dynamical restrictions. 

(B) Since time in Newtonian cosmology t is absolute, all clocks of any observers in 
an inhomogeneous Newtonian universe are always synchronized. So arranging the same 
measurement time intervals T by such a covering is undoubtedly possible. The situation 
with a space region S G with a volume Vs determined by the measurement device 
resolution is more delicate. A typical argument against a possibility to arrange a covering 
with space region S of the same shape and volume makes use of the fact that an evolving 
cosmological fluid configuration may expand or shrink during its evolution. However, one 
can easily show that this argument does not even appear if all physical hypotheses have 
been clearly stated and checked. Indeed, in accordance with the Newtonian cosmological 
macroscopic hypothesis (|[) the characteristic scale V^^^^ of the space regions S must be 

— 1/3 

always much less than a characteristic cosmological scale L5, Vg << Ls- Therefore 
any peculiarities of the cosmological fluid evolution of a characteristic cosmological scale 
Ls cannot affect the measurement device resolution and change either the shape or the 
volume of the averaging space regions S placed at different positions at different times. 
Therefore, such a covering by the same space regions S is always possible. 

(C) The values of the space-time averaged fluid field {f{x\t))grp calculated for such a 
covering depend now on T and Vs as free parameters and such an average field becomes 
effectively a local single- valued function of the reference time t and position {x^} only, 
and the correlation function {f{x\t)h{y^ , s))grp becomes a local single- valued function of 
each pair of variables. This fundamental result and its consequences for the analytical 
properties of the space-time averages (j^) are analyzed in the next Section. 

6 The Reynolds Conditions for Average Fluid Fields 

In choosing some particular averaging schemes, in addition to the conditions for the con- 
struction of the averages fluid fields, see Section |, one must also investigate a possibility 
to formulate the general requirements which such a scheme should have. From the point 
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of view of Newtonian cosmology, the most important of these general requirements is, 
of course, that the application of an averaging rule to the Navier-Stokes-Poisson equa- 
tions will allow one to get sufficiently simple equations for the average cosmological fluid 
fields. Reynolds has found such a set of conditions while he was using time averages for 



the derivation of the averaged Navier-Stokes equations [16 . These conditions have been 



reformulated later for any averaging procedure applicable in fluid mechanics |T^, |^ 



p8| , in the classical macroscopic electrodynamics |2l[]-[^ and in general relativity |JT2 



r5[ , P5| . Let us consider fluid quantities / = f{x\t) and h = h{x\t) which characterize 
an inhomogeneous and/or anisotropic cosmological fluid conflguration evolving in time 
such that both fluid flelds have nonvanishing temporal and spatial derivatives (|^). The 
Reynolds conditions are formulated |1^, [Q, p8| as follows. 



The Reynolds Conditions For any cosmological fluid fields f{x^,t) and h{x\t) the 
space-time averages {f{x\t))gj, and {h{x\t))gj, must satisfy the following conditions: 

(i) the space-time averaging is a linear operation 

{af{x\t) + bh{x\t))g^ = a{f{x\t))g^ + b{h{x\t))g^, if a, 6 = const, (16) 

(a) the space-time averaging commutes with the partial differentiation 

^^{f{x\t))sT = {p{x\t))sT, ^{f(^\t))sT = {^J{x\t))sT, (17) 

(Hi) the idempotency of the space-time averages 



[{f{x\t))sTh{y\s))^^ = {f{x\t))sT{h{y\s))^^ 



or 



\{f{x\t))sT)sT=^f(^'^t))sT. 



One can show that the space-time average (|^) satisfles the Reynolds conditions, see 
121 , |1T3| , [p5|] for discussion and proofs of relevant issues. 



Theorem 1 (The Reynolds conditions for the space-time averages) The space- 
time average ^ is a linear averaging procedure satisfying the Reynolds condition (^7^. // 
a covering ( [7^ j hy the averaging time intervals T and the space regions S is determined 
through the cosmological fluid configuration, then 

(1*) the space-time average field {f{x\t))^rp is a local single-valued function of the 
reference time t and the position {x^}, which depends on the value ofT and Vs as param- 
eters, such as, 

^ «)>.,.^0. (^^-4^V/(.'.0>,,^0. (19) 



dx^ dt dtdx^J ^ ' ' / ST ' \dx^ dx^ dx^ dx' J ' / st 
(2*) the space-time average field {f{x'^,t))gj, satisfies the Reynolds condition (jl^. 
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Proof. By its definition the space-time averaging procedure is a linear operation 
(plGl). Let us now assume that in accordance with the First and Second conditions of 
Section |^ a covering (|15D has been arranged for every moment of time t at every position 
{x*}. Then one can determine an averaged field {f{x^,t))grp and define now a spatial 
partial derivative of the average field along the coordinate axis x*^*^ as 

^ {f(^'^^))sT = J?Lo + ^^'''')>^T - (/(^^^)>.J (20) 

where the average value (/(x* + Ax'-*-', t))^^ is taken at the reference point {x*' + Ax'^'^-'} 
and at the moment of time t. Using the definition of the space-time average ex- 
panding /(x* + Ax^*-*, t ) into a Taylor series around the {x*^} and taking into account the 
conditions ([T5|) on the averaging time intervals T and space regions S, one easily derives 
the commutation formula for the averaging and the spatial derivative, 

^(/(^^^))^T = (^/(x^^)).^, (21) 

which does not contain explicitly either T or Vs as a consequence of the covering (|T^). 
Defining a time derivative of the averaged field {f{x\t )) gj, in the same way as (]20|), one 
can derive the commutation formula for the averaging and the time derivative, 

|(/(x^^))5T = (|/(x^t))5T, (22) 

which does not contain explicitly either T or as a consequence of the covering (p!5|). 
Calculation of the second derivatives of the averaged field (/(x*,t))^^ in accordance with 



the definition (^) bring the fundamental result The formula means locality of 

an average field {f{x\t))grp, that is, it is a local single-valued function of the reference 
time t and position {x^}, which does depend on the value of T and Vs as parameters. One 
can work with these averages as with usual functions depending on time t and position 
{x'^} by differentiating and integrating them, expanding into series, solving differential 
and integral equations for them using standard techniques of mathematical physics. The 
Reynolds condition (|T^ is given by (|21|) and (H). QED 

A particular choice of the measurement (averaging) time interval T and the mea- 
surement resolution (averaging) space region S can be fixed at any convenient stage of 
analysis. All formalism and the form of the averaged equations do not depend on the 
choice of values T and Vs- When this choice is made, the average field {f{x\t))gj, ac- 
quires particular values at each point {x'^} at each instant of time t, which refiects the 
physical conditions of the measurements and the cosmological fiuid dynamics. 

The last Reynolds condition (^) is known to be an approximate condition |T^, |n 



p5| which holds in Newtonian cosmology due to the Newtonian cosmological macroscopic 



hypothesis (H). 

Corollary 1 (The idempotency of the space-time averages) The space-time aver- 
ages for a cosmological fluid configuration satisfying the Newtonian cosmological macro- 
scopic hypothesis are idempotent and the Reynolds condition ( |7^ j holds. 
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It should be pointed out here that the physical meaning of the idempotency ( ]TBD of the 
space-time averages (^) is that under the Newtonian cosmological macroscopic hypothesis 
(^) the measurement time interval T and a space region S G with a volume 



s 



determined by the measurement device resolution guarantee that the errors made to hold 
the conditions (|T8D are always less than the errors of the time and space measurements, 
0{Xt/T) and 0{Xs/Vg ^''^) (see [0 and references therein). 



A satisfactory formal analysis of the idempotency property ( [I8| ) of the space-time, 
time and space averages is still lacking, though it is extensively used for the derivation 
of the averaged equations and analysis of various fluid distributions. It has been shown, 
however, that there are certain classes of volume averaging procedures which satisfy all 
Reynolds conditions (p!^)-(|T8|), but these procedures and corresponding averaging kernels 
do not have a clear physical interpretation and they did not find direct applications in 
fluid mechanics (see [p^], pOl for a discussion and references). 



7 The Properties and Conditions for the Correlation 
Functions 

One can now establish the analytic properties of the correlation function {f{x\ t)h{y\ s))gj, 
(|ll]) and find its properties analogous to the Reynolds conditions (p!^)-([T8|). 



Corollary 2 (The properties of the space-time correlation functions) If a cover- 
ing ^Tdj) by the averaging time intervals T and the space regions S is determined through 
the cosmological fluid configuration, then the two-point second order moment function 
(fH|j, {f{x\t)h{y^,s))gj,, of two cosmological fluid fields f{x\t) and h{x\t) has the fol- 
lowing properties: 

(1*) it is a bilocal single-valued function of the reference times and positions, {x^,t) 
and {y^,s), which depends on the value of T and Vs as parameters, such as all second 
antisymmetrized derivatives with respect to all pairs of the variables {x'^,t]y^,s) vanish, 

(2*) the two-point second order moment function {f{x\t)h{y^ , s))gj, satisfies the con- 
ditions of the partial differentiation, 

d ^ , ^\ /dfix\t) 



dx^ 



f{x\ t)hiy\ s))^^ = {^%^hiy\ .)^^^ , 

(/(x^t)%^.)>^^ = (|^^^^%^.)^^^, etc (23) 



ST 

„fc +. „,k 



for all variables (x , t; y , s), 

(3*) the two-point second order moment function {f{x\t)h{y\s))grp is idempotent 

{{f{x\t)h{y\s))^^g{z\r))^^ = {f{x\t)h{y^,s))^^{g{z\r))^^ or (24) 

{{fix\ t)h{y^, s))^^)^^ = {f{x\ t)h{y^, s))^^ , (25) 
where g = g{z'^, r) is another cosmological fluid field. 
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Therefore one can work with the space-time correlation functions as with usual local 
functions of many variables depending on reference times and positions by differentiating 
and integrating them, expanding into series, solving differential and integral equations for 
them using standard techniques of mathematical physics. 

It should be noted here that the conditions (|23| ) mean that a partial differentiation 
of a correlation function produces a correlation function of the same order including the 
corresponding partial derivative of a fluid field. 

There is an important asymptotic property 
Cgj,(x\t;y^ , s) (PI^D, which is of particular interest in Newtonian cosmology. 

The asymptotic condition for the central two-point second order space-time 
correlation function For any two cosmological fluid fields f{x^,t) and h{x^,t) which are 
correlated ( |7^ and satisfy the Newtonian cosmological macroscopic hypothesis their 
central two-point second order moment function C gj^{x^ ^t]y\ s) tends to zero when the 
distance between the points {x*} and and the time interval between t and s infinitely 
grow 

Cfl{x\ t; y^, s)^0 as [6ij{x' - y'){x^ - y^)] ^ oo, \s-t\^oo. (26) 



of the central correlation function 



This property has a very clear physical meaning stating that cosmological fluid fields 
are almost independent at extremely remote points in space and for remote moments of 
time. So the central correlation functions ([13D for the cosmological fluid fields in the New- 
tonian universes of Class LC, see Section ^ always satisfy the asymptotic condition (p^). 
It should be pointed out, however, that for the Newtonian universes of Classes L and C 
(1^) and (§) the central correlation functions ( [T^ ) may require different asymptotic condi- 
tions. The asymptotic conditions on the correlation functions (|lT]) should be considered 
each time separately for each Newtonian universe. 



8 The Time Averages of a Cosmological Fluid 



Let us consider now the time averages JT^, , which are a particular case of the 
space-time average (^. A time average is defined as a limit of successive measurements 
of a fluid field / = f{x'^,t) over the measurement time interval T until the fluctuations in 
its average values become acceptably small. 

The time average of a fluid fleld / 



t2 



f{x\t + t')dt'. 



(27) 



Here the time average (^) is a well-defined function of a reference moment of the abso- 
lute time t to which the time average {f{x^,t))j, is prescribed, t & T, and the position 
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{x*} where the measurements of f{x^,t) are carried out by an observer with a measure- 
ment device resolution region S, {x*} G S. The time average (|27|) satisfies the Reynolds 
conditions ([l|)-(|lD. 

Corollary 3 (The Reynolds conditions for the time averages) // a covering ([7^ 
by the averaging time intervals T and the space regions S is determined through the cos- 
mological fluid configuration, then the time average (/(x*,t))^ is a local single-valued 
function ( |7^ j of the reference time t and the position {x^} depending on the value of T 
as a parameter and satisfying the Reynolds conditions ((7^ j-/f7^j; 
(1*) the time averaging is a linear operation 

{af{x\t) + bh{x\t))^ = a{f{x\t))^ + b{h{x\t))^, if a, 6 = const, (28) 

(2*) the time averaging commutes with the partial differentiation 

= (|/(x\t))T, = {^fi^\t))T, (29) 

(3*) the time average is idempotent 

{{f{x\ t))Th{y\ s))^ = {f{x\ t))T {h{y\ s))^ or {{f{x\ t))^)^ = {f{x\ t))^. (30) 

Time averages can be applied for the Newtonian universes when a cosmological fiuid 
distribution satisfies the Newtonian cosmological macroscopic hypothesis (^) for the tem- 
poral variation scales (A^, Lt) and there is either a local fiuid field spatial fiuctuation scale 
A5, or a characteristic cosmological spatial scale L5. In the first case the cosmological 
fiuid configuration is completely inhomogeneous in space with a typical spatial variation 
scale A5 of local fiuid field fiuctuations, to have the vanishing space-time average (H) of a 
fluid fleld f{x\t) over a space region S C with a volume Vs such as A^ << Vg ^^'^ , 

{f{x\t))^^ = 0, {f{x\t))^=^0. (31) 

In the second case the cosmological fiuid configuration is completely spatially homoge- 
neous ( ^91 ) on the cosmological spatial scale L5. Then the space-time average of a 
cosmological fiuid field f{x\t) over a space region S C with a volume Vs such as 
Vg « Ls does not change it (p!8D to smooth f{x\t) effectively only over the time 
interval T, 

{fix\t))^^ = {fit))^. (32) 
Let us now consider the condition when the time average for a cosmological fluid flelds 



/ = /(x*,t) does not depend on time. The following important theorem takes place |r7 |. 



Theorem 2 (The time average over an infinite time interval) For any bounded 
cosmological fluid quantities f{x\t) and h{x\t) the time average ^Ij ) over an infinitely 
large time interval T does not depend on time 

lunjf{x\t))^ = {f{x'))^. (33) 
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The two-point second order time correlation function {f{x\t)h{y^, s))rp defined due to ( [7i[ j 
as 

{f(x\t)h{y\s))^ = ^ j f{x\t + t')h{y\s + t')dt' (34) 

ti 

and evaluated over an infinitely large time interval, depends on the difference s — t only 

^lim {f{x\ t)h{y\ s))^ = B{x\ y^- s-t). (35) 

Proof. The proof is as follows. Let us consider the same time average over a measure- 
ment time interval T at the reference point t, {f{x\t))rp, and at another reference point 
s, {f{x\s))rp, s,t & T, and estimate the difference {f{x\s))j, — {f{x\t))j,. Due to the 
idempotency (|30|) of the times averages this difference is approximately zero since under 
the Newtonian cosmological macroscopic hypothesis (^ the measurement time interval 
T guarantee that the errors made to hold the conditions (|30| ) are always less than the 
errors of the time measurements, 0{\t/T), see Section ^ However, now one would like 
to find out a condition when the idempotency holds precisely for the time average defined 
as (p7l). Assuming s > t calculation of the difference {f{x\s))rp — {f{x\t))rp gives 



{f{x\s))^-{f{x\t)) ^ 



'T 



T 



S+t2 t+t2 

f{x^,r)dr~ j f{x^,r)dr 

\s+t^ t+tl 



(36) 



For a bounded function of time f{x^,t) the difference ( p^ ) becomes infinitely small as 

I T 1^ oo, 

{fix\t))^ = {f{x'))^ as T-. 00, (37) 

which proves (^). The proof of (|35|) is straightforward is one considers the limit of the 
correlation function {f{x\t)h{y\ s))^ as T — > oo 



Bix\ y^; s-t)= ^lim {f{x\ t)h{y^ , s))^ = 

t2 

lim ^ [ f{x\ t + t')h{y^, {s-t)+t + t')dt' (38) 

T^oo T J 
h 

which shows that {f{x\t)h{y^,s))j, averaged over an infinitely large time interval can 
only depend on (s — t), but not on s and t individually. QED 

The condition (^) can be also easily understood if one observes that a time average 
over the whole interval of time allowed by the evolution of a Newtonian universe must 
not depend on time. Indeed, an imaginary measurement of a time dependent physical 
quantity over the whole time domain where the quantity is determined gives a number 
provided such a time average value is well-defined. The time correlation function (|3^ can 
depend in this case only on the time difference between the correlated values of the fluid 



19 



fields ( P^D since the average fields themselves does not depend on time (P^D- The higher- 
order time correlation functions can be shown all depend only on the time differences, for 
example, B{x\ , z^; s — t,r — t) = limT^oo {f{x\ t)h{y^ , s)g{z^, '^))t' 

A fluid configuration satisfying everywhere the following conditions for averaged fluid 
fields and their time correlation functions over a finite measurement time interval T 

(/(x\t)>^ = (/(xO>^, {fix\t)hiy\s))^ = Bix\y^;s-t), etc (39) 

is called stationary. 

The physical meaning of stationarity is transparent since all the equations governing 
such a fluid motion are time- independent. 

The Newtonian stationary universes A cosmological fluid configuration is stationary 
if the fluid fields satisfy the conditions at all positions {x''} during its evolution. 

It has been shown in Section II-8 that a particular class of the stationary Newtonian 
universes, namely, a static homogeneous and isotropic Newtonian universe does not exist 
if the Newtonian cosmological constant vanishes, A = 0. 



9 The Space Averages of a Cosmological Fluid 



Let us consider now the space averages jl^ which are another particular case of the space- 



time averages @. A space average is defined as a measurement of a fluid field / = f{x\t) 
over a space region S G with a volume Vs determined by the measurement device 
resolution such as the fluctuations in its average value become acceptably small. 

The space average of a fluid field / 

{fi^\t))^ = ^^ [ f{x' + x'\t)dV' (40) 



s 



Here the space average (^0]) is a well-defined function of a reference position {x^} to which 
the space average (/(x*,t))^ is prescribed, {x^} G S, and the moment of time t when the 
measurements of /(x*, t) are carried out by an observer with a measurement time interval 
T, t E T. The space average (HOl) satisfies the Reynolds conditions (p!6|)-([T8|). 



Corollary 4 (The Reynolds conditions for the space averages) If a covering i^T3{ ) 
by the averaging time intervals T and the space regions S is determined, then the space 
average (/(x*,t))^ is a local single-valued function (jl^ of the reference time t and position 
{x^} depending on the value ofVs as a parameter and satisfying the Reynolds conditions 

m-m-- 

(1*) the space averaging is a linear operation 

{af{x\t)+hh{x\t))g = a{f{x\t))^ + h{h{x\t))^, if a, 6 = const, (41) 
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(2*) the space averaging commutes with the partial differentiation 



(3*) the space average is idempotent 

{{f{x\t))sh{y\s))^={f{x\t))s{h{y\s))^ or {{f{x\t))s),= {f{x\t))s. (43) 

Space averages can be applied for the Newtonian universes when a cosmological fluid 
distribution satisfies the Newtonian cosmological macroscopic hypothesis (|^) for the spa- 
tial variation scales {Xs,Ls) and there is either a local fluid field temporal fluctuation 
scale At, or a characteristic cosmological temporal scale Ly. In the first case the cosmo- 
logical fluid configuration is completely inhomogeneous in time with a typical temporal 
variation scale Xt of local fluid field fluctuations to have the vanishing space-time average 
(H) of a fluid field f{x\t) with a time measurement interval T such as << T 

{f{x\t))^^ = 0, {f{x\t))^^0. (44) 

In the second case the cosmological fluid configuration is stationary (|39|) on the cosmolog- 
ical spatial scale Lt- Then the space-time average (^) of a cosmological fluid field f{x\t) 
over a time measurement interval T such as T << L5 does not change it (pIS]) to smooth 
f{x\t) effectively only over a space region S 

{f{x\t))^^ = {f{x^))^. (45) 

Let us now consider the condition when the space average for a cosmological fluid 
field / = f{x\t) does not depend on the reference point {x*} . The following important 



theorem takes place |T7 



Theorem 3 (The space averages over the whole space) For any bounded cosmo- 
logical fluid quantities f{x^,t) and h{x^,t) the space average ^4^) over the whole space 
S = does not depend on a position {a;*} 

hm {f{x\t))^ = {f{t)),. (46) 

The two-point second order space correlation function {f{x^,t)h{y^,s))g defined due to 
([Hi as 

{f[x\t)h{y\s))^ = ^ [ f{x' + x'\t)h{y^ + x'\s)dV' (47) 

Js 

and evaluated over the whole space S = depends on the difference y^ — only 

hm {f{x\t)h{y\s))^ = B{y^ -x';t,s). (48) 
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Proof. The proof can be easily accomplished in a completely analogous manner as the 
proof of the similar properties for the time averages (p3D and (B^). QED 



The condition (^61) can be also easily understood if one observes that a space average 
over the whole space allowed by the definition of a Newtonian universe must not 
depend on spatial coordinates. Indeed, an imaginary measurement of a position dependent 
physical quantity over the whole space where the quantity is determined gives a number 
provided such a space average value is well-defined. The space correlation function { ^7\) 
can depend in this case only on a difference between the correlated values of the fluid 
fields (^ ) since the average flelds themselves do not depend on a space position (^). 
The higher-order space correlation functions can be shown all depend only on the time 
differences, for example, B{y^ — x\ — x^; t, s, r) = limvg^oo {f{x\ t)h{y^ , s)g{z^, r))^. 

A fluid conflguration satisfying everywhere the conditions 

{f{x\t))^ = {f{t)),, {f{x\t)h{y\s))^ = B{f-x'-t,s), etc (49) 

for the average fluid characteristic flelds and their space correlation functions over a space 
measurement device resolution region S G is called homogeneous. The physical 
meaning of homogeneity is transparent since all the equations governing such a fluid 
motion are independent of fluid positions {x*}. 

The Newtonian homogeneous universes A cosmological fluid configuration is ho- 
mogeneous if the fluid fields satisfy the conditions at all positions {x'^} during its 
evolution. 

It has been shown in Section II-8 that a particular class of homogeneous Newtonian 
universes, namely, a static homogeneous and isotropic Newtonian universe does not exist 
if the cosmological constant vanishes, A = 0. The class of homogeneous and isotropic 
Newtonian cosmologies has been studied in Sections II-6 and 11-12. 

It is important to compare the formal deflnition (^91) with the formulation of the New- 
tonian cosmological principle, see Section II-6, which deflnes the homogeneity condition 
for a cosmological fluid conflguration in terms of the data measured by a privileged family 
of observers. The homogeneity conditions (^) set precise constraints on the character of 
the averaged fluid flelds as they are measured by the observers endowed with measure- 
ment devices having the same measurement time interval T and the same measurement 
device resolution determined by a space region S (|T5|). 



10 The Space-time Averaging out of Material Deriva- 
tives 

In derivation of the averaged systems of the Navier- Stokes- Poisson equations (II-4), (II-8), 
(II-9) and (11-10), or the Navier- Stokes-Poisson equations in terms of kinematic quantities 
(II-4), (II-8), (11-12), (II-51)-(II-56), (II-59)-(II-61), (I-64)-(II-66), (11-69), (11-28) and (II- 
70) one needs to consider carefully the problem of taking an average value of the Navier- 
Stokes fleld operator, that is, the left-hand side of Eq. (11-10). As it has been pointed out 
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in the Section ^ an averaging out of this operator (|I|) does not give the material derivative 
of an averaged fluid velocity. Let us consider now the space-time averages (u\x^ ,t)^ 
(^) of the cosmological fluid velocity u'^x^.t) and make use of the Reynolds conditions 
(|16D-(|18D after taking the space-time average of the Navier-Stokes field operator, that is, 
the material derivative of the fluid velocity. One arrives at the expression 



ST 

at ^^^^i^ (50) 

which again is not equal to the material derivative of the averaged velocity t))^^ 
because of the presence of the single-point second order space-time correlation function 
of the fluid velocity and its spatial derivative (^u^ {x^ ,t)uj{x'' ,t)'^ Therefore, a formula 
for performing averaging out the Navier-Stokes field operator is definitely very important 
as it will affect the structure of the field operator of the averaged Navier-Stokes equation 
and finally that of the averaged Navier-Stokes- Poisson equations. 

In an inhomogeneous Newtonian universe any fluid particle has different values of its 
instant velocity u'^{x\t), the density p{x\t) and the pressure p{x^,t) as it moves along 
its path. A cosmological fluid configuration with the space-time average fields for the 
velocity {u\x^ ,t)) gj,, the density {p{x\t))gj, and the pressure {p{x\t))grp may have a 
different dynamics as compared with the original fluid configuration since the average 
fields are not equal to the original fiuid fields (H) in general, see Section |[ Indeed, a fluid 
particle moving in an inhomogeneous self-gravitating fluid conflguration with the velocity 
u^{x\t) has an equation of motion (I-l), 

x' = x\e,t), (51) 

as a solution to the initial value problem (^. Due to the Newtonian cosmological macro- 
scopic hypothesis (H), an averaging region x T is considered effectively as a single 
"point" for the averaged cosmological fluid fleld {f{x\t))grp, see Section ^. Therefore, 
one can consider now the same fluid particle as moving in the averaged self-gravitating 
fluid cosmological conflguration with the average velocity {u^i^x^ ,t)) ^j, such as (|^) 

{u\x\t))^^^u\x\t). (52) 

The equation of motion 

X'=X\i\t) (53) 

of the fluid particle moving in the averaged cosmological fluid conflguration is different 
from the equation of motion (|5T|) for the fluid particle moving in an inhomogeneous 
cosmological conflguration with the velocity u\x'',t) because a solution to the initial 
value problem (^), 

-^ = {u\x\t))^^, X\0)=e, (54) 
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is different from tliat to (|^) for tlie same initial values x^{0) = ^* and X*(0) = So 
the paths of the same fluid particle as it is considered moving within an inhomogeneous 
cosmological fluid distribution and within the corresponding averaged cosmological fluid 
distribution have to be different in general. 

A fluid quantity / in the framework of the field description is a function of the Eulerian 
variables, / = /(x*, t), and it is also a function of the material variables, / = /(^*, t), such 
as 

f = f[x'{e,t),t] or f = f[eix\t),t], (55) 

where functions 

e = e(x^t) (56) 

define the initial position {C}, C = x*(0), of the fluid particle which is at any position 
{x*} at a moment of time t, see Section 1-5. The average fluid field (/) = {f{x\t))grp is 
now a function of the Eulerian variables (1531) and the material variables 



e = e{x',t), (57) 

such as 

{f) = {f)[X\e,t),t] or {f) = {f)[eiX\t),tl (58) 



the functions (^) being related by Eqs. (|53D and (|57D considered as the laws of change 



of variables. Then the change in the quantity (/) in course of the averaged fluid motion 
can be characterized by the two different time derivatives 

at at |X^=const at Ot |e=const 

The partial derivative d {f) jdt gives the rate of change of (/) with respect to a fixed 
position {X*}, while the material, or convective, derivative d (/) jdt measures the rate of 
change of (/) with respect to a moving fluid particle. 

The material derivative of the position of a fluid particle defined by the second Eq. 
( |59D for (/) = X* is called its average fluid velocity (m*) 

(«^> (e,t) = ^ = . ' («^> (^^'^) = («^> (^''^) = («^(^^^)>5T• (60) 

dl Oh |g'=const 



The notion of the average fluid velocity (m*) ( |60D plays a fundamental role in the field 
formulation of the fluid mechanics of the averaged flows in the same way as the fluid 
velocity vJ" plays a fundamental role in the field formulation of fluid mechanics of the 
inhomogeneous flows, see Section II-5. 

With the definitions of time derivatives (|5^ one can calculate the material derivative 
of any averaged quantity (/) as 



dt 



|(/>(e.o 



|{/) |A"(e,«),(| 

|$»=const ) |5'=const 



dt |X^=const dX^ 



dX\i\t) 
dt 



(61) 

|{*=const 
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that can be written in the following form by using the average fluid velocity defined by 
Eq. (|0D: 

m = m^u)m or m=m+u)m (62) 

dt dt ^ \ / dX^ ° dt dt + / dx^ ■ ^ ^ 

This formula relates the material and spatial derivatives through the average velocity 
field in the field picture of the average fluid motion and it expresses the rate of change 
in (/) with respect to a moving fluid particle located at a position {X*} at time t. It is 
completely analogous to the material derivative of the fluid field f{x\ t) (1-12), see Section 
1-6, 

df df ,df 

-r = ^ + UTr-, (63) 
dt dt dx' ^ ^ 

where it relates the material and spatial derivatives through the velocity field u'^{x\t) in 
the field description of fluid motion, see Section 1-5, and it expresses the rate of change 
in f{x^, t) with respect to a moving fluid particle located at a position {a;'} at time t. 

The Eulerian coordinates (X*,t) of the averaged cosmological fluid configuration cor- 
respond to the Eulerian coordinates (x*, t) of an inhomogeneous cosmological fluid config- 
uration since the averaged value {f{x\ t))grp has been prescribed to a reference time t and 
a position {x*} which are the same moment of time t and the same position {X*} = {x*}. 
One must distinguish between the equations of motion of fluid particles moving in an inho- 
mogeneous fluid distribution and in the corresponding average fluid distribution averaging 
out particular relations, but one can always use the spatial coordinates {x*} instead of 
{X*} in the final averaged formulae (|5^). 

One can now prove the following important theorem establishing the rule for space- 
time averaging of the material derivatives. 

Theorem 4 (The space-time average of the material derivative) // a covering 
j^Idj ) by the averaging time intervals T and the space regions S is determined through the 
cosmological fluid configuration, then the space-time average of the material derivative 
of a cosmological fluid field f{x^,t) is given by the formula 

dt/„ \m dx'/„ dt ^\ dx'/ ^ ' dx' ' 



where 



is the single-point second order space-time correlation function of the fluid velocity u^{x\ t) 
and its spatial derivative df{x\t)/dx\ 

Proof. The proof is straightforward by taking the space-time average (^ of the material 
derivative of a cosmological fluid field / = f{x\t) with using the Reynolds conditions (p!7|), 
the definitions of the average cosmological fluid field (/) (|58D and the average cosmological 
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fluid velocity (m*), the definition and expression for the material derivative of the average 
fluid field (|60|) and ( |62D and the introduction of the single-point second order space-time 
correlation function (^). QED 

There are two equivalent useful forms of the formula 



df\ d{f) / ,df 



d{f) 
dt 



cm 



u 



dl 



u 



(67) 



where [u\ is the central single-point second order moment function of the fluid 
velocity u'^{x^,t) and the spatial derivative df{x\t)/dx^ of the fluid field f{x\t). 

The formulae (|6^)-(|67|) are of fundamental significance in the derivation of the averaged 
Navier-Stokes-Poisson equations. An immediate result is their application to the space- 
time averaging out of the field operator of the Navier-Stokes equation (??). 

Corollary 5 (The averaged Navier-Stokes field operator) If a covering ^idj) by the 
averaging time intervals T and the space regions S is determined through the cosmological 
fluid configuration, then the space-time average (j^ of the field operator of the Navier- 
Stokes equation (11-10) has the following form: 



du" 
'dt 



ST 



dt ^ ^ dx^ 



ST 



d{u') 
dt 



+ u 



dx^ 



dx^ 



where 



dxi 



U-' 



dxi 



(.^t)^(.^(x^t)M^) 



(68) 



(69) 



ST 



is the single-point second order moment function of the fluid velocity u''{x^,t) and its 
spatial derivative du''{x^,t)/dx^ . 

Thus, the averaged Navier-Stokes field operator has changed its structure under the 
space-time averaging (|^) as it has been suggested by a heuristic argument (|50D . The 
averaging rule ( |68D plays a central role in establishing the form of the averaged Navier- 
Stokes-Poisson equations. 



Averaged Material Derivatives and the Reynolds 
Transport Theorem 



11 



The derivation of the formulae (|6^ ) - (|67|) for a space-time averaging of the material deriva- 
tives of cosmological fluid fields has required a thorough analysis of the averaged fluid 



motion, see Section 10. The foundation of these formulae lies in the possibility to arrange 
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the covering (]T3|). Then in accordance with the First and Second conditions, see Section 
^ for a covering of Newtonian space-time by the averaging time intervals T and the space 
regions S typical in some defined sense, one can determine the space-time average field 
{f{x\t))gj. which is a local single-valued function of the reference time t and position 
{x'^}, which depends on the value of T and Vs as parameters and satisfies the Reynolds 
conditions (|T^)-(0 ), see Sections ^ and ^. The covering (^) has been shown to be 
consistent with the Newtonian cosmological macroscopic hypothesis @ and it associates 
a standard measurement device with the measurement time interval T and the measure- 
ment resolution space region S with each point {x^} at each instant of time t. This is the 
physical foundation of the Reynolds conditions which are always assumed for derivation 
of the averaged equations and relations for space-time, time, space and ensemble averages 
in fluid mechanics Jl^, ||20l, [^, the Maxwell macroscopic electrodynamics |2l|]-[^ and 



general relativity [Q, ||T3|, ^ 

Let us now consider a covering of Newtonian space-time by the measurement time in- 
tervals T and the space regions S G determined by the measurement device resolution. 



which are taken to be typical in different than the covering (|T^) sense. If one requires a 
more general covering 

f-' 

at each moment of time t and position {x''} that means T = const and Vs 7^ const. 
In fact, it is actually impossible to require a covering other than T = const with the 
two first conditions (ffOD since the time is absolute in Newtonian cosmology. A choice 
of space regions S G such as Vs 7^ const with the two last conditions (|70D held 
is possible formally. Such a choice associates a measurement device with its resolution 
space region S at each point {x''} at each instant of time t with its characteristic scale 

—1/3 

Vs changing in dependence of a position in space. In order to have Vs 7^ const, the 
Newtonian cosmological macroscopic hypothesis should be violated by two possible 
ways, namely, either A5 ~ V5. << L^, or A5 << — Ls- In the first case the 

change in the space regions S is due to a spatial scale A5 of the local cosmological fluid 
field fluctuations. It cannot be taken seriously from the physical point of view since 
observers would be affected by the peculiarities of the fast and sharply changing local 
fields and any space averaging cannot be performed. In the second case the change in the 
space regions S is due to a spatial cosmological scale Ls and the observers can perform 
space measurements in principle. However, the measurement space resolution scale Vg 
is of the same order as a spatial mean cosmological scale that implies, as a matter of fact, 
that the space averaging does not resolve any cosmological scale at all. Thus, a covering 
( |70D should result in violation of the Newtonian cosmological macroscopic hypothesis 
(P) and the impossibility to determine a family of observes with physically reasonable 
measurement devices which are able to perform proper space-time measurements of the 
cosmological fluid characteristics. 

The thorough analysis of the averaged fluid motion, physical foundations underlying 
the definition of averages and the formalism for establishing their properties, see Sections 
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^ ^ and |TU], has resulted in the derivation of the formulae (|61D-(|^ for the space-time 
averaging out of the material derivatives of cosmological fluid fields. It is especially im- 
portant for Newtonian cosmology where all these issues have a direct relation to the way 
how our Universe is observed, how the cosmological data are measured and how a theoret- 
ical model takes into account the measurement device characteristics and the conditions 
of experiments. Using a covering different than ( ]T3| ) and/or using the definitions of a 
space-time, time or space averages different than (||), ( ^Tf ) and ( ^OD would lead to a vio- 
lation of the Newtonian cosmological macroscopic hypothesis and an arrangement of 
a physically unreasonable set of measurement devices. That would affect, as a result, the 
analytic properties of the averages to abandon the Reynolds conditions (|l6l)-(|l^) and the 



formulae (|6^) 



It is therefore important to analyze possible consequences of such different choices in 
order to make a comparison with the approach presented here for the formulation of the 
foundations of Newtonian cosmology and the averaged Navier-Stokes-Poisson equations. 
Let us consider here a possibility to use the Reynolds transport theorem (1-20) for deriva- 
tion of the formula for a space averaging out of the material derivatives. On the basis of 
the Reynolds transport theorem is easy to show the following: 

Corollary 6 (The modified Reynolds transport theorem) The rate of change of 
the ratio of the volume integral F(t) (1-19), 



F{t)= J fix\t)dV, (71) 

m 

of the cosmological fluid field f{x\t) over an arbitrary closed fluid region S(t) moving 
with the fluid to the value Vj^it) of the volume of the region (1-24), 

if) (^) = ^ = T7^ / fdy^ (72) 



is given by 



Vj,{t) V,,{t) 

m 



''^"'-fflwW/l^lw-mwIl^wo^ (73) 



dt \dt j \ dx^ j [ dx'^ 

Proof. The formula evidently follows from the Reynolds transport theorem (1-20) and 
a corollary of the theorem (1-25) giving the rate of change of the volume (1-24) of the 
region E(t). QED 

The function of time {/} (t) measures the ratio (|7^ ) of the value of the quantity 



/(x*,t) inside an arbitrary closed fluid region S(t) moving with the fluid to the value 
of the region's volume. The modified Reynolds transport theorem (^) determines the 
change of the ratio {/} (t) rather than of the integral F{t) ( |7TD as the Reynolds transport 
theorem (1-20) does and they are completely equivalent. All results in analysis of the fluid 
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kinematics made on the basis of (1-20) can be achieved by using the modified Reynolds 
transport theorem ([f3|). 

Though the quantity {/} (t) ( [72D looks similar to a space average {f{x\t))g (^OD of 
the cosmological fiuid field f{x^,t), it cannot be taken as a proper definition of a space 
average of f{x\t) due to the following reasons: 

(#1) the quantity {/} (t), as well as the quantity F{t), are well-defined functions of t 
only, since they are defined only for arbitrary closed fiuid regions moving with the 
fiuid by the meaning of the Reynolds transport theorem (1-20); 

(#2) the quantity {/} (t) cannot serve as a space average, since no arrangements have 
been made for the construction of a fiuid field {/} (t) determined at each point {x*} and 
each moment of time t, and, therefore, for an arrangement of a set of measurement devices 
throughout the cosmological fiuid configuration, see the First and Second conditions of 
Section ^ 

(#3) the regions S(t) moving with the fiuid change their shape and volume, which 
raises all the problems with ( [7D| ) discussed above; 

(#4) the Reynolds transport theorem (1-20) is known to be an integral statement fully 
equivalent to the Euler expansion formula (1-15), see [jT| and references therein, therefore 
any conclusions gained on the basis of application of the modified Reynolds transport 
theorem ([73|) are completely equivalent to the corresponding statements for the integrand 
of (O). 

0^5) the Reynolds transport theorem (1-20) and the modified Reynolds theorem ( [73|) 
do not take into account by their physical meaning and their conditions that the fiuid 
particles of an averaged cosmological fiuid configuration have different equations of motion 
X* = X\^^ , t) ([53|), see Section [ly, as compared with the equations of motion = x*(^'' , t) 
(|5TD of the fiuid particles of the corresponding inhomogeneous fiuid configuration and the 
region S(t) is supposed to move along the fiuid particles' paths determined by ([^); 

The modified Reynolds theorem ([73|) cannot be therefore used as a formula for a space 
averaging out of the material derivatives instead of the formulae (p^)-(p7|). 



12 The Ensemble Averages of a Cosmological Fluid 



The use of space-time, time or space averages (^), (^) and (|40| ) is very convenient from 
the practical point of view, but leads sometimes to analytical difficulties as it has been 
discussed above. It is desirable therefore to find another method of defining the average 
fiuid fields such that it would have well-defined analytic properties and would be universal. 
A convenient definition of this type is known to be provided by the statistical picture of 



the fiuid fields as random fields (see 17 , [IS 



which has been established by Richardson 



and reference therein), 
35[| , Taylor |^, Kolmogorov ||3^ 



and Kampe de Feriet p9 



The basic feature of the statistical approach to the motion of fiuid is the transition 
from the consideration of a single moving fiuid configuration given by the equations of 
motion of fiuid particles x* = x^{C,-',t) ([5l| ) or ^"^ = C,^{x^,t) (p6D to the consideration of 



29 



the statistical ensemble of all similar moving fluid configurations created by some set of 
fixed initial and boundary conditions. If the values of the fluid velocity field /(x*, t) for an 
evolving inhomogeneous fluid configuration are measured many times under the same set 
of conditions, the measured values (/(x*, t)) for the same moment of time t and the same 
position {x*} are expected to be different due to local spatial and temporal fluctuations of 
/(x*, t). However, the arithmetic mean of all these values at the moment of time t and the 
position {x*} can be expected to tend to a definite value if a sufficiently large number of 
measurements has been carried out. This arithmetic average over all possible realizations 
of the fluid quantity f{x\t), that is, all the values {f{x\t)) as measured, in principle, by 
an observer, is called an ensemble average of f{x\t) at the time t and the position {x*}. 

The ensemble average, or the probability mean, of a fluid field / 

+00 

{f)Ei^'^t)= J fP^,.,){f)df. (74) 



Here the probability density function P(x\t){f) defined for all times t and points {x*} of 
the cosmological fiuid configuration determines the probabilities of different values /(x*, t) 
among all measured as 



+00 



^(xSt)(/)c?/ = Probabity[/</(x\t)</ + rf/], J P(,.,,)(/)rf/ = 1. (75) 

~oo 

The quantity Probabity [/ < f{x\ t) < f + df] in (^) stands for the probability of f{x\ t) 
to have values between / and f + df . A fiuid field / having a definite probability density 
is called the random variable and the set of all possible probabilities P(x\i){f' 1 f") = 
Probabity [/' < f\x\t) < /"] corresponding to / is called its probability distribution. 

The ensemble average (0) is a well-defined function of a reference moment t of the 
absolute time, t and a reference position {x^}, {x^} e S, where the ensemble average 
(/)_£; (^N^) is evaluated by a measurement device with the measurement time interval T 
and the space region S G determined by the measurement device resolution. It should 
be pointed out here that the measured values of / may be time, space, or space-time 
averages in dependence on the physical conditions of a cosmological fiuid configuration 
under consideration. The nature of measurements must be taken into account when 
making physical predictions and interpretations on the basis of the statistical treatment 
of fiuid dynamics. 

Because the Navier-Stokes-Poisson equations are nonlinear one faces necessity to deal 
with averaging products of the fiuid fields. The ensemble average (|74[) for a product 
of cosmological fiuid fields / = /(x*,t) and h = h{x'^,t) defines the so-called ensemble 
correlation function. 

The two-point second order ensemble correlation, or moment, function of fluid 
fields / and h The two-point second order correlation function is the ensemble average 
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of a product of f{x^,t) and h{y\ s) 



{fh)^{x\t;y\s)= fhP^^.,t)iy.,s)if,h)dfdh. (76) 



Here P(x\t)(yi ,s){f ih) is the two-dimensional joint probability density function defined 
for all times t and s and points {a;*} and {?/*} of the cosmological fluid configuration 
determines the probabilities of different values f{x\t) and h{y^,s) among all measured 
as 

Pix^,t)(y^,s){f, h) = Probabity [/ < f{x\ t) < f + df, h< h{y^ , s) < h + dh], 

+00 

P{x\t){y^,s)if, h)dfdh = 1. (77) 

The quantity Probabity [/ < f{x\ t) < f + df, h < h{y^ , s) < h + dh] in (|77D stands for 
the probability of f{x\ t) and h{ii\ t) to have simultaneously values between / and f + df 
and between h and h + dh, correspondingly. 

The correlation function (^) is a well-defined function of the reference moments t and s 
of the absolute time, t,s E T, and of the reference positions {x''} and {y''}, {x''}, {y^} G S, 
with the ensemble average (fh)^ {x\ t; y\ s) being prescribed to the pair of the Eulerian 
variables (x',t) and {y-^,s) and symmetric in these variables, 

{fh)j,{x\t;y^,s)E={hf)j,{y^,s;x\t). (78) 

From the physical point of view the two-point ensemble average ( |76D corresponds 
to a simultaneous measurement of / and h by an observer with a measurement device 
with the time measurement interval T and the measurement device resolution region 
S G E^. As far as the Newtonian cosmological macroscopic hypothesis (^ is satisfied for 
a cosmological fiuid configuration under consideration a choice of the reference times t 
and s and the reference positions {x*^} and {y'^} is arbitrary within T and S. 

The correlation functions of fiuid fields are fundamental characteristic functions re- 
sponsible for essentially nonlinear phenomenon in evolving fiuids, such as, for example, 
turbulence and hydrodynamic instability [Q, [1^, ^Q\, [^]. If some of such functions do 



not vanish for a fiuid configuration, then the corresponding fiuid fields are said to have 
correlations. 

The fluid field correlations A cosmological fluid configuration has a central two-point 
second order ensemble correlation, or moment, function {x^,t;y\ s) , if there are, at 
least, two cosmological fluid fields f{x^,t) and h{x\s) such that 

CP {x\ t; y\ s) = ifh)^ {x\ t; y\ s) - (/)^ {x\ t) {h)^ {y^ , s) ^ 0, (79) 
C'i\x\t;y^,s) = CP{y\s;x\t) (80) 
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in some open space region U C E for an interval of time At, x ,y E U, t,s E At. 

One can define the central multi-point higher-order ensemble correlation functions, 
C^^\x^,t;y\ s; z'^ ,r) and so on, by using higher-dimensional joint probability density 
functions P{xi,t){yi ,s)(z3 ,r){f, h, g) and so on, whenever it is necessary for analysis of the 



dynamics of a Newtonian universe by making definitions similar to (^) and (79 ). 

In order to characterize the ensemble averaged fluid configuration one needs to de- 
termine now the ensemble averaged fields {f)^{x\t) and the correlation function fields 
(fh)^ (x*, t] , s) for each moments of time t and s at all positions {x^} and {y^}. From 
the physical point of view, that means that it is possible, at least in principle, to carry out 
the measurements of the quantities / and h for all instants of time during the evolution 
of a fluid configuration and at its each space position. This requires a covering by the 
averaging time intervals T and the averaging space regions S to be determined throughout 
the cosmological fluid configuration in accordance with the First and Second conditions, 
see Section |[ It will be assumed that a covering (p!5D has been arranged throughout the 
cosmological fluid configuration and the Newtonian cosmological macroscopic hypothesis 
(H) is valid for this set of measurement devices. 

Thus one can formulate the main hypothesis on the statistical description of Newtonian 
universes. 

The Hypothesis of the statistical nature of Newtonian universes: The cosmo- 
logical fluid fields of the components of the fluid velocity vector u'^{x\t), the fluid density 
p(a;^*,t), the fluid pressure p{x\t) and the Newtonian gravitational potential (f){x\t) are 
assumed to be random fields defined for all times t and positions {x*} of the cosmolog- 
ical fluid configuration with the corresponding probability densities P(-j.j P(^^i i){p), 
P{x\t){p) O'^d P(^x\t){4>)- If the cosmological fluid fields have correlations (fT^j, that is, 
they are statistically interconnected, it is assumed that there exist the joint probability 
densities for corresponding fluid fields. 



13 The Properties of the Ensemble Averages 

The ensemble average (|74D satisfies the Reynolds conditions (p!^)-([T^). 

Corollary 7 (The Reynolds conditions for the ensemble averages) The ensemble 
average (/)^ {x\t) is a local single-valued function 

f d d d d \ f d d d d \ 

of the reference time t and position {x^} and it satisfies the Reynolds conditions ([7^^- (17^.- 
(1*) the ensemble averaging is a linear operation 

{af + bh)j,{x\t) = a{f)j,{x\t) + b{f)j,{x\t), if a, 6 = const, (82) 
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(2*) the ensemble averaging commutes with the partial differentiation 
(3*) the ensemble average is idempotent 

[{f)^{x\t)h)^{y\s) = {f)^{x\t){h)^{y\s) or ((/))^ (x\ t) = t). (84) 



Proof. The formulae (|8TD-(p4D follow immediately from the definition of the ensemble 
average QED 

The analytic properties of the ensemble averages are the same as those of space-time, 
time and space averages. It is remarkable fact that they all are rigorous properties for the 
ensemble averages. This is one of the fundamental advantages of using this type of mean 
fluid fields. 

The formula for the ensemble averaging out of the material derivatives can be proved 
as follows. 

Theorem 5 (The ensemble average of the material derivative) The ensemble av- 
erage i fm ) of the material derivative of a cosmological fluid field f{x'^,t) is given by the 
formula 

dt/^~\dt^'' dx^/^~ dt +\ 9x7^ ^""'E dx^ ^^^^ 



where {u^f^i)^ = {u^f^i)^ {^''^t) is the single-point second order correlation function ^/dj ) 
of the fluid velocity u''{x^,t) and the spatial derivative df{x\t)/dx\ and the material 
derivative of (f) ^ {x\t) is 

dt dt ^ ' dx' ^ ' 

Proof. The proof is straightforward by taking the ensemble average ( |7^ of the material 
derivative (|63D of a cosmological fluid fleld / = f{x\t) with using the Reynolds conditions 
(|83D , the deflnitions of the ensemble average {f)^{x\t) (^) of the cosmological fluid 
fleld and the ensemble average {u^)^ {x^, t) of the cosmological fluid velocity u^i^x^ , t), and 
the introduction of the single-point second order correlation function {u^f.i)E (^'')^) (^)- 
QED 

There are two equivalent useful forms of the formula ) which are completely anal- 
ogous to the formulae (|66|) and 



Corollary 8 (The ensemble averaged Navier-Stokes field operator) The ensem- 
ble average j^T^ ) of the field operator of the Navier-Stokes equation (11-10) has the following 
form: 



dt \dt dx^/ ^ dt \ dx^ / ^ ^ dx^ 
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(87) 



where (u^u^^^^ = (u^u^^^^^x^^t) is the single-point second order correlation function 



( fTqj of the fluid velocity u'^{x\ t) and its spatial derivative du^{x , t)/dx^ and the material 



derivative of the ensemble averaged fluid velocity {u"^) ^ {x\t) is given by 

One can now establish the analytic properties of the ensemble correlation function 
{fh) ^ {x\ t] y^ , s) (^) and prove the properties analogous to the conditions (|23D -(|25 D . 



Corollary 9 (The properties of the ensemble correlation functions) The two-po- 
int second order moment function ^Tdj ), (fh)^ {x\ t; y^ , s), of two cosmological fluid fields 
f{x\t) and h{x\t) has the following properties: 

(1*) it is a bilocal single-valued function of the reference times and positions, {x^,t) 
and {y^, s), such as all second antisymmetrized derivatives with respect to all pairs of the 
variables {x^, t; y^, s) vanish; 

(2*) the two-point second order moment function {fh)^{x'^,t;y^,s) satisfies the con- 
ditions of the partial differentiation 

^ ifh)^ {x\ t- y\ s) = (^^)^ t; y\ s), 

^^{fh)j,{x\t-y\s) = (^^h^ {x\t-y\s), etc (88) 



for all variables {x^, t; y'', s); 

(3*) the two-point second order ensemble moment function (fh)^ [x^, t; y\ s) is idem- 
potent 

{ {fh)^ {x\ t; y\ s)g)j^ (^^ r) = {fh)^ {x\ t; y\ s) {g) ^ {z\ r) or (89) 

{{fh))j,{x\t;y^,s) = {fh)j,{x\t;y\s) (90) 
where g = g{z'',r) is another cosmological fluid field. 



Proof. The formulae (p8D-(pOD follow immediately from the definition of the ensemble 
correlation function (^). QED 

Thus, one can work with the ensemble correlation functions as with usual local func- 
tions of many variables depending on the reference times and positions by differentiating 
and integrating them, expanding into series, solving differential and integral equations for 
them using standard techniques of mathematical physics. 

It should be noted here that the conditions (|88| ) mean that a partial derivative of an 
ensemble correlation function produces an ensemble correlation function of the same order 
including the corresponding partial derivative of a fluid field. 

There is an important asymptotic property of the central ensemble correlation 



(2') 

function {x^,t; y^ , s) (^), which is of particular interest in Newtonian cosmology. 

The asymptotic condition for the central two-point second order ensemble 
correlation function For any two cosmological fluid fields f{x'',t) and h{x'',t) which 
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are correlated ( [7^ and satisfy the Newtonian cosmological macroscopic hypothesis 
their central two-point second order moment function C^\x^,t;y^ , s) tends to zero when 
the distance between the points {x*} and {y^} and the time interval between t and s 
infinitely grow 

dE\x\t;y^,s) ^0 as [6ij{x' - y'){x^ - y^)^^^ ^ oo, |s-t|^oo. (91) 



14 The Ergodic Hypothesis of Newtonian Cosmology 

The ensemble averages are undoubtedly more convenient for analytic description of mov- 
ing fluids. Adopting a hypothesis on the existence of the probability distributions for 
all fluid fields, such the Hypothesis of the statistical nature of Newtonian universes in 



Newtonian cosmology, see Section one can use the mathematical techniques of the 
modern probability theory and statistical mechanics. The ensemble averages ( |7iD and 
(|76D are defined uniquely (0) throughout a cosmological fluid configuration and have all 
properties naturally required, such as the Reynolds conditions (]82|)-(p^ and the condi- 
tions (^)- (pO|) . It is important, however, to note that in the framework of this approach, 
an important new question arises regarding the comparison of the theoretical deductions 
following from the statistical treatment of the dynamics of fiuids with the data of direct 
space-time, time and space measurements. 

According to its definition (^) the ensemble average is a mean value taken over all 
possible values of the fiuid field under study, which are measured by a measurement device 
with a time measurement interval T and a measurement space region 5* determined by 
the device resolution. Thus, to determine the average values by an experiment or, say, by 
cosmological observations, with an acceptable high accuracy, one should carry out a large 
number of measurements in many set of repeated similar observations or experiments. In 
practice, however, no such multiply repeated measurements are available, or even possible 
to perform. Therefore in most cases one is forced to determine the average values from 
the data taken in the course of a single experiment where the data obtained either by 
a space-time, time or space averaging over a time interval T and/or a space region S. 
From the physical point of view it means that the assumption about the existence of 
the probability distributions, see the Hypothesis of the statistical nature of Newtonian 
universes, see Section O, does not by itself eliminate the problem of the validity of 



using the space-time (H), time ( p?]) or space averages (PUD for description and analysis of 
moving inhomogeneous fluids. But in this context the problem has acquired a different 
formulation and interpretation. Indeed, instead of investigating the particular properties 
of these averaging procedures, one should establish how close the space-time, time or space 
averaged values of the fluid flelds measured in real observations and experiments are to the 
corresponding ensemble average values. It is known from the statistical mechanics that a 
replacement of an ensemble average value of a random quantity calculated over all possible 
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states of this quantity by the directly measurable time average of the same quantity may 
be only possible if the averaging time interval T becomes infinitely large. Under this 



condition the time average converges to the corresponding ensemble average ^8 



pO| , [p^ , p^ , p3| . In certain very special cases, the validity of this assumption may 
be proved rigorously. In the most cases of interest this fundamental physical hypothesis 
of the consistency between theoretical models, their experimental verification and the 
real physical phenomenon remains unproved and it is usually adopted as an additional 
hypothesis called the Ergodic hypothesis. It is formulated here for Newtonian cosmology. 

The Ergodic hypothesis of Newtonian cosmology: // a cosmological fluid configu- 
ration satisfies the Newtonian cosmological macroscopic hypothesis ^ and the Hypothesis 
of the statistical nature of Newtonian universes, see Section and a covering ( |73| j by 
the averaging time intervals T and the space regions S is determined throughout the cos- 
mological fluid configuration, then for any cosmological fluid field f{x^,t) the following 
conditions are assumed to hold: 

(1) its ensemble average (/)^ {x'',t) ( f^D is equal to its time average {f{x\t))rp over 
an infinitely large time interval T and both averages do not depend on time t (|55D, 

lnnjf{x\t))^ = (/(xO>^ = (92) 



(2) its ensemble average (/)^ {x\t) jj^ ) is equal to its space average {f{x\t))g over 
the whole space S = and both averages do not depend on a position {x*} (^), 



hm {f{x\t))^ = {f{t)),= {f)^{t), (93) 

— >oo 



(3) its ensemble average {f)^{x\t) ( |7^ is equal to its space-time average {f{x\t))gj, 
(j£) over an infinitely large time interval T and over the whole space S = E^ and both 
averages do not depend on on time t and a position {x*} ^3^ ) and (U^), 

[f{x\t))^^={f),^ = {f)^. (94) 



lim lim 



Thus, if proved, the Ergodic hypothesis of Newtonian cosmology (p2|), (p^) and ( p4D 
guarantees that the values of cosmological fluid quantities measured by an observer in 
a Newtonian universe by means of space-time, time and/or space averagings are equal 
within the accuracy of measurements to the values predicted theoretically on the basis of 
the statistical treatment of the evolution of the Newtonian universe. 



15 The Averaged Navier-Stokes-Poisson Equations 
in Kinematic Quantities 

The system of the averaged Navier-Stokes-Poisson equations in kinematic quantities can 
be now established. Since the space-time (^ and ensemble averaging, (H) and (0), have 
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the same analytic properties ( [T^ ) and ([ST|), satisfy the Reynolds conditions (|T7p-(0) and 



!2| )-(|8^), and the same formulae for averaging out the material derivatives (0), (|68D and 
(|85D , (pTl), the notation (/) will used hereafter for both averaging procedures. Both 

procedures lead to the same form of the averaged equations, though they are physically 
equivalent, strictly speaking, only under the Ergodic hypothesis of Newtonian cosmology 
(P^, see Section |H[ Use of the time and space averages (P7D, (^), see Sections ^ and 
can be made every time when the physical conditions of a Newtonian universe require 
this and the corresponding parts of the Ergodic hypothesis of Newtonian cosmology, (B^p 
and (p3D, are valid. 



Theorem 6 (The averaged Navier-Stokes-Poisson equations) // a covering ([71 
by the averaging time intervals T and the space regions S is determined through the 
cosmological fluid configuration, then the averaged system of the Navier-Stokes-Poisson 
equations in terms of the kinematic quantities (11-4), (^^'6), (H-S), (11-12), (II-51)-(II- 
56), (II-59)-(II-61), (1-64)- (11-66), (11-69), (11-28) and (11-70) is as follows. 

The averaged Raychaudhuri evolution equation for the averaged expansion 
scalar {6) 



dt 



'i"'^)- W-^l + 2«-^> - H-4^G (p) -AS" (A.),, = 0, (95) 



The averaged propagation equation for the averaged shear tensor (aij) 



1 



5.,(2 {a') + {u') - 5^' {A,)^,) + {uj^uj,) + (E.,) - {A^^ = 0, (96) 



The averaged propagation equation for the averaged vorticity vector (u;*) 

d{u;') 



dt 



+ 



The averaged kinematic decomposition of the averaged tensor {uij) 
(ui) J = (cxij) + ^6,j (6) + (uij) , {(Xij) = (aji) , (6) = 6'^ {a,j) , {uij) = - (uji) , (98) 



The averaged first identity for the averaged tensor {uij) 
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The averaged second identity for the averaged tensor {uij) 

i^i) ,jk - i^i) ,kj , (100) 

The averaged first constraint equation 

^^■'((^..),.-K),.)-^(^),=o, (101) 

The averaged second constraint equation 

{^'), = 0, (102) 

The averaged third constraint equation 

(5^(^((a,,)^, + (a;,,)_,)£-)^^ = 0, (103) 

The averaged first integrability condition 

<bO, = ^(p),<. (104) 

The averaged second integrability condition 

sOfc^ ^O, (105) 



The averaged third integrability condition 



The averaged equation of continuity as the evolution equation for the averaged 
fluid density (p) 
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The averaged equation of state 

(P) = (pip)) , 



(108) 



The averaged total acceleration (Ai) by the averaged Navier-Stokes equation 

{A,) = - (^-^p.'j , (109) 

The averaged Poisson equation for the averaged Newtonian gravitational po- 
tential (0) 

6'^ {g,) = AnG (p) - A, or 5^^' (0) = 47rG (p) - A, (110) 



The averaged identity for the averaged Newtonian gravitational acceleration 

{9^) 

^^^■'(^.). = 0, (111) 



The averaged tidal force tensor (Eij) for the averaged Newtonian gravitational 
potential (0) 



The averaged evolution equations for the averaged tidal force tensor {Eij) 
d{Eij) / f^ dEij \ /„,fc\ d {Ejj) skris,im I TP \ 

35"-' {<yk(:,Ej)i) - 5^=' (^^(.S,-)^) + 47rG (pa,,) = 0, (113) 



Proof. The averaging out of the system of the Navier-Stokes-Poisson equations in terms 
of kinematic quantities (II-4), (II-6), (II-8), (11-12), (II-51)-(II-56), (II-59)-(II-61), (1-64)- 
(11-66), (11-69), (11-28) and (11-70) is straightforward by applying the Reynold conditions 
(P^-(P^ or (p2D-(p4D, and the formulae for averaging out the material derivatives (0), 
(§§ or (|S3|), (|7p. QED 

The system of the averaged Navier-Stokes-Poisson (^ )-( |113| ) contains the following 
independent set of the single-point second order moment functions: 
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The single-point second order moments in the averaged Navier-Stokes-Poisson 
equations: 

The fluid velocity with the spatial derivatives of the kinematic quantities aij, 9 and 
uj'^, the fluid density p and the tidal force tensor Eij 




(114) 



The kinematic quantities a^j, 9 and uj'^ themselves 

{9'), 6'' , {9a,,), (uiCUj) , {9lj') , 5^' {a.^uj^) , (115) 

The fluid density p with the kinematic quantities aij and 9 and the pressure gradient p^i 



{p9) , {paij) , 




The tidal force tensor with the kinematic quantities aij, 9 and 

{9Eij) , 5ij5^'^5^"' {akiEran) , S''^ {o-k{iEj)i) , 5^^ {uk{iEj)i) . (117) 



The system of the averaged Navier-Stokes-Poisson (^ )-( |113D is not closed with respect 
to these single-point second order moment functions. The system should be supplemented 
by the equations for these moment functions ( |114| )-( [TT7| ). A derivation of a system of the 
equations for the second moments is known 



16- 1 



26 



3, F 



to produce 

the higher-order moments of fluid fields. This is the so-called closure problem. To make 
an analysis possible and also take into account properties of particular fluid conflgurations 
one must formulate some conditions for the fluid fleld moments of some order to terminate 
this inflnite set of equations. 
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